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We study the effect of fluids on the angle of repose and the segregation of granular matter poured 
into a silo. The experiments are conducted in two regimes where: (i) the volume fraction of the 
fluid (liquid) is small and it forms liquid bridges between particles thus giving rise to cohesive forces, 
and (ii) the particles are completely immersed in the fluid. The data is obtained by imaging the 
pile formed inside a quasi-two dimensional silo through the transparent glass side walls and using 
color-coded particles. In the first series of experiments, the angle of repose is observed to increase 
sharply with the volume fraction of the fluid and then saturates at a value that depends on the size 
of the particles. We systematically study the effect of viscosity by using water-glycerol mixtures to 
vary it over at least three orders of magnitude while keeping the surface tension almost constant. 
Besides surface tension, the viscosity of the fluid is observed to have an effect on the angle of repose 
and the extent of segregation. In case of bidisperse particles, segregation is observed to decrease and 
finally saturate depending on the size ratio of the particles and the viscosity of the fluid. The sharp 
initial change and the subsequent saturation in the extent of segregation and angle of repose occurs 
over similar volume fraction of the fluid. Preferential clumping of small particles causes layering to 
occur when the size of the clumps of small particles exceeds the size of large particles. We calculate 
the azimuthal correlation function of particle density inside the pile to characterize the extent of 
layering. In the second series of experiments, particles are poured into a container filled with a fluid. 
Although the angle of repose is observed to be unchanged, segregation is observed to decrease with 
an increase in the viscosity of the fluid. The viscosity at which segregation decreases to zero depends 
on the size ratio of the particles. 

PACS number(s): 45.70.Mg, 47.50.+d, 64.75.+g, 81.05.Rm 



I. INTRODUCTION 

The presence of small amounts of liquid can have a 
considerable influence on the properties of granular mat- 
ter. For example, the angle of repose of a wet granular 
pile is greater than that of a dry pile made of the same 
material. Cohesive forces are introduced because of liq- 
uid bridges that are formed between particles. Cohesivity 
causes jamming in the flow of granular matter even in a 
regime where dry granular matter may flow. Wet gran- 
ular matter is also observed to segregate less than dry 
granular matter as particles cannot move easily relative 
to each other. Although these qualitative facts are well 
known, a detailed knowledge of the properties even in 
comparison with dry granular matter has not yet been 
attained. 

The properties of dry granular matter have attracted 
a considerable number of studies. Granular matter in 
a box which is tilted does not flow until a certain an- 
gle is exceeded at the surface. This angle usually called 
the maximum angle of stability 9 m depends on the fric- 
tional properties and the packing of the granular matter. 
However, if the granular matter is in motion, it can re- 
lax down to a lower angle which is called the angle of 
repose 9 r . The existence of the lower 9 r may be due to 
the lower kinetic friction compared to static friction be- 
tween the particles and other dynamical effects [Q. The 
observed hysteresis may also occur because there are in- 



finitely many meta-stable states above 9 r in which the 
system can get trapped . 

In case of a very large pile of granular matter, it has 
been argued that 9 m is unaffected by the cohesivity intro- 
duced by the presence of the fluid. This is true provided 
that the frictional properties between the particles are 
unaltered due to the presence of a lubrication layer. The 
basic argument is that while frictional force which helps 
hold a pile together increases with the size of the pile, 
the cohesive force is constant and becomes insignificant 
in comparison However the effects can be impor- 
tant for a laboratory or industrial scale system. Recently, 
there have been a number of studies that have tried to re- 
late the fluid content and its surface tension to observed 
increase in angles. 9 r measured by Schiffer and collabo- 
rators using the draining crater method is observed 
to initially increase linearly and then fluctuate around 
a constant value. For an ideal sphere-sphere contact, in- 
creasing the liquid content will lead to decreased cohesive 
force and hence a lowered 9 r . Therefore they proposed 
that particle surface roughness plays an important role 
and a cone-plane type of contact between particles may 
be more appropriate. They also believed that the fluid 
coats the particles evenly. Mason et al. || have studied 
the distribution of the fluid on the surface of the parti- 
cles using optical fluorescence microscopy and scanning 
electron microscopy. These studies appear to show that 
the fluid is in fact trapped in the menisci that are formed 
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at the asperity on the grain surface. They have therefore 
claimed that the behavior of 9 m at low volume fractions 
of the fluid is consistent with a scaling, theory based on 
the surface roughness of the grains MM. Other studies 
note that considerable differences between 9 r and 9 m can 
arise when humidity is present and have shown that the 
wetting properties of the fluid can be also important |j) . 

The identification of the occurance and the extent of 
segregation is of considerable importance in material han- 
dling in industry. Preferential percolation of small parti- 
cles through layers of larger particles is usually identified 
as leading to segregation when granular flows occur at a 
free surface [|lO],[ll) . Stratification can occur in such sys- 
tems if rough and smooth shaped grains are present . 
The pre sence of the fluid alters the percolation of parti- 
cles |Q and thus has a profound influence on the progress 
of segregation in granular matter. It is thus surprising 
that the effects of the fluid on the progress of segrega- 
tion has not been investigated in detail till recently. In 
a recent publication, we reported that the presence of a 
fluid drastically reduces segregation in bidisperse granu- 
lar matter that is poured into a silo E3. 

Although it is clear that the surface tension of the fluid 
is important in determining the cohesive force and hence 
6 r , the role of the viscosity of the fluid v may be less ob- 
vious. It is a common observation that a finger dipped in 
honey feels sticky in comparison to a finger dipped in wa- 
ter. This occurs because stronger forces are required to 
move the liquid into regions that are vacated. Therefore 
viscous forces may have a significant effect on the dy- 
namics of the grains. It is thus important to investigate 
the effect of v on both 9 r and on size separation. Fur- 
thermore, viscous effects grow in importance for smaller 
dimensions. It may be possible to understand the ef- 
fects of humidity on powders, that are smaller in size in 
comparison, by investigating effects of viscosity on segre- 
gation in granular systems. This may be useful because 
quantitative studies with powders are more difficult to 
conduct. 

In this context, we report a detailed study of the an- 
gle of repose and segregation of granular matter in the 
presence of a fluid using high resolution digital imaging. 
Some of the results have been reported in a Letter 
and here we report additional data and analysis. We first 
investigate the effect of the size of the particles on 9 r as a 
function of the volume fraction of the fluid. We also dis- 
cuss in detail our experimental observation of the effect 
of the viscosity v of the fluid on 9 r . We then study the 
extent of segregation by visualizing the color-coded glass 
particles using bidisperse glass beads with size ratio r. 
Interestingly, the decrease in segregation is observed to 
saturate at similar values of volume fraction as 9 r . It can 
be also noted that small particles preferentially clump in 
a bidisperse mixture and this leads to subtle effects on the 
nature of the particle spatial-distribution as the volume 
fraction of the fluid is increased. As the volume fraction 
of the fluid increases and as the limit is approached where 
all the interstitial volume is saturated with the fluid, one 



can expect that the effects of the cohesion introduced 
by liquid bridges to become less important. This may 
have an effect on 9 r and size separation. However, there 
are practical difficulties in investigating high volume frac- 
tions because the fluid tends to drain. Therefore we also 
conducted experiments when the particles are completely 
immersed in a fluid to understand this limit. 



II. EXPERIMENTAL APPARATUS 

A schematic of the experimental apparatus is shown in 
Fig. [lj A rectangular silo of dimensions 50.0 cm x 30.5 cm 
and a width w of 3.0 cm is used for the experiments. The 
flow is visualized through the glass side-walls of the silo 
using a 1000 pixel x 1000 pixel Kodak ES 1.0 digital 
camera. The glass particles and fluids used are listed in 
Table | and Table FJ respectively. The granular sample 
is prepared in batches by throughly mixing 1.0 kilogram 
of each of the two kinds of particles with the fluid before 
the mixture is placed inside the reservoir. The volume 
fraction Vf of the fluid is calculated as the ratio of the 
volume of the fluid to the volume of all the particles. The 
volume of the particles corresponds to the weight of the 
particles divided by the density of the glass beads which 
is 2.4 gem -3 . 

The wet granular material is first filled into a reservoir 
and then drained through a pipe into the silo. In pre- 
vious experiments we demonstrated that segregation can 
occur during draining of a wide silo [^5| , and therefore 
a tall cylindrical reservoir is used to minimize segrega- 
tion from occuring during pouring. Pipes with various 
diameters are used to control the flow rate Q. Larger 
diameters are used at higher Vt in order to maintain 
Q ~ 2.2 g s _1 . The reservoir is raised at a slow con- 
stant rate with a stepper motor and a system of pulleys. 
The slow upward velocity of the reservoir allows the par- 
ticles to accumulate inside the pipe before flowing down 
the surface and reduces the kinetic energy of particles 
due to free fall. Such precautions are necessary other- 
wise particles acquire substantial kinetic energy during 
free fall and tend to bounce back from the surface sev- 
eral times leading to size separation and stratification in 
the case of bidisperse particles. Smaller particles bounce 
back higher from the surface of the pile and land further 
down the pile in contrast to what is usually observed, 
i.e. small particles at the top of the sandpile. This ef- 
fect coupled with periodic avalanching observed at the 
surface can lead to alternate layers of mixed and small 
particles. By using the pipe, we ensure that the inter- 
action of particles is restricted to the surface of the pile. 
Thus the number of mechanisms involved in segregation 
is reduced and the system is simpler. 

We varied the time that the sample was mixed and 
found it had little effect on 9 r . We also waited a vary- 
ing amount of time before pouring to test the effect of 
evaporation. We found that only in the case of water, 
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there was a small influence if we waited for one hour for 
volume fractions below Vf < 2 x 10~ 3 . Therefore, to min- 
imize the effects of evaporation, the mixture was poured 
immediately after mixing. 



III. THE EFFECT OF INTERSTITIAL FLUID ON 
THE ANGLE OF REPOSE 

A triangular shaped granular pile is formed after the 
material is poured into the silo. The flow of the granu- 
lar material is continuous for dry particles, but becomes 
increasingly stick-slip as Vf is increased. The flowing re- 
gion is roughly 10 layers deep for dry particles, but the 
flow region becomes deeper and is not as well defined as 
Vf is increased. The resulting pile for monodisperse glass 
particles is shown in Fig. || and it can be observed that 
the surface of the wet pile is at a greater inclination than 
the dry pile. 

A. Size dependence of the angle of repose 

9 r is measured by obtaining images of the sand pile 
after the silo is filled and fitting a straight line to the 
surface of the pile. The obtained slope is averaged over 
10-15 images as a centimeter deep layer of granular mat- 
ter is added and a new surface is created. 9 r is plotted 
as a function of Vf for particles with various diameters d 
in Fig. ||(a). 9 r is first observed to increase sharply with 
Vf and then saturate approximately. The value of Vf at 
which 9 r begins to saturate is V c . The saturation value 
of 9 r is observed to increase with the size of the particles. 
This fact can be qualitatively noted from Fig. |j. 

The increase and approximate saturation of 9 r is quali- 
tatively similar to those earlier observed using the drain- 
ing crater method @J{|. Quantitative differences exist 
perhaps because of the difference in the geometries. As 
discussed in the introduction, both the increase in 9 r and 
saturation is surprising if one assumes ideal spheres in 
contact. To explain this property, work has focused on 
the fact that particles are not perfectly smooth but have 
asperities. Using Mohr-Coloumb analysis and making as- 
sumptions for the height and distance over which the sur- 
face of the sphere fluctuates, Halsey and Levin Q| found 
that 9 m can increase linearly and then saturate to a value 
given by: 

tanfr m — sec tan (k)\ (1) 

dp g gH 

where, k is coefficient of static friction, T is the surface 
tension of the fluid, g is acceleration due to gravity, p g 
is the density of the particles, and H is height of the 
pile. It has been assumed here that the capillary force 
F ca p between the particles is given by irTd. They as- 
sume that the saturation will occur approximately when 
the wetting is determined by the macroscopic curvature 



of the particles, which is greater than the scale of the 
roughness of the surface of the particles. According to 
the model, volume V c is given by \ 2 d where A is length 
scale over which fluctuations occur. The model predicts 
that 9 m decreases with d, which is consistent with our 
observations. Although there is a qualitative agreement 
between our data for 9 r and the model for 9 m , there is 
considerable disagreement related to the value of V c . For 
example for 1 mm particles, where the particle surface 
fluctuations are about 1/im, the model predicts satura- 
tion for V c ~ 10~ 5 . However the observed V c is signifi- 
cantly higher and is approximately 5 x 10~ 3 . 

While it is possible that the disagreement can arise 
because m is calculated and 9 r is measured, additional 
effects may have to be taken into account in calculat- 
ing the average cohesive force due to capillary bridges. 
For example, it has been assumed that only the capillary 
forces between particles in contact is important. How- 
ever it is possible that as Vf increases, particles that 
are slightly away from each other can also form liquid 
bridges. Therefore, it is possible that the average num- 
ber of liquid bridges per particle increases with Vf . This 
effect may be important in determining the regime over 
which 9 r can increase. 

We plot in Fig. |^(b), the 9 r for particles with d — 0.5 
mm and d — 3.1 mm and the bidisperse mixture (BD-3) 
to illustrate the behavior of 9 r when particles of different 
diameters are present. The 9 r for the bidisperse mix- 
tures are located in between the 9 r of their components. 
A similar effect is observed for the other size ratios. 

Besides the increase of 9 r due to the addition of the 
fluid, the surface of the pile is observed to become rough 
[see Fig. ||(b),(d)]. In addition, the roughness of the sur- 
face of the pile is greater for the smaller particles. We 
plot the mean square of the deviation of the surface from 
a straight line \ 2 f° r the particles mixed with water for a 
fixed Vf — 24 x 10~ 3 in Fig. ||(c) to quantitatively show 
the increase in roughness with decreases in particle size. 
We observe that smaller particles clump more than larger 
particles at similar Vf . We therefore believe that the in- 
creased roughness of the surface is because of clumping 
of the particles. As the particles form more clumps, the 
surface becomes increasingly rough and \ 2 grows. The 
standard deviation of 9 r obtained by repeating measure- 
ments for a particular Vf was reported in Ref. || as a 
measure of clumping. The reported observations are con- 
sistent with our more direct measurement of the surface 
heights. We will later see that the preferential clumping 
of small particles has an important effect on the progress 
of segregation and introduces layering. 



B. Effect of viscosity of the fluid on the angle of 
repose. 

To investigate the effect of viscosity on the flow and 
the angle of repose, we kept the surface tension of the 
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fluid constant and change v of the fluid. Because water 
and glycerol have very similar surface tension, we used 
mixtures of these two fluids to change v |16|. Glycerol 
dissolves in water and the mixture is homogenous after 
it is stirred for a few minutes. Using water-glycerol mix- 
tures we obtained 9 r as a function of Vf . 

In Fig. ||(a) wc plot 6 r for d = 0.5 mm beads as a 
function of Vf, for water, glycerol, and a water-glycerol 
mixture. Fig. ^(a) shows that 6 r increases with v. The 
increase is observed to be sharper and the saturation is 
observed to occur at higher values for glycerol compared 
to water. We chose Vf = 24 X 1CP 3 in the saturated 
regime to plot 9 r as a function of v for monodisperse par- 
ticles and bidisperse mixtures in Fig. ||(b). We observe 
that saturation value for 9 r depends on v and increases 
with v. 

We have measured 9 r up to 48 hours after the granular 
mixture is poured into the silo. While a few local rear- 
rangements occur up to a few minutes after pouring, we 
do not find measurable changes in the surface, thus in- 
dicating the absence of creep. We have thus shown that 
the viscosity of the fluid clearly influences the observed 
8 r . In the next section we will discuss the viscous force 
between the particles due to the fluid. 



spheres. They showed that Eq.(4) fits their experimental 
data. 

The value of the viscous force in our experiments can 
be estimated assuming that ^f- is of the order of the ve- 
locity of the particles flowing down the surface. For 1 mm 
particles, if ho is of the order of 10 /im then F V i s ~ 10~ 7 
N for water and F vis ~ 2 x 10~ 4 N for glycerol. The 
capillary force estimated using F cap — nTd for 1 mm 
particles is 10~ 4 N. Therefore the viscous force for glyc- 
erol is relevant whereas for water the viscous force may 
be insignificant. These improved results are also consis- 
tent with our previous estimates obtained neglecting the 
curvature of the particles . 

Because viscous force decreases with relative velocity 
it appears surprising that the viscosity of fluid plays any 
role at all in determining 9 r in our experiments and we 
return to this point later in the paper. However it is 
clear that such forces are important in determining the 
extent of segregation as it occurs when the particles are 
in motion. 



IV. EFFECT OF INTERSTITIAL FLUID ON SIZE 
SEGREGATION 



C. Estimates of the viscous force 

Let us consider two particles that are moving away 
from each other. A schematic diagram is shown in 
Fig. |(a). Fig. |(b) shows 0.1 mm particles mixed with 
glycerol for Vf = 24 x 10~ 3 . Based on Reynolds lubrica- 
tion theory, an expression for the viscous force has been 
derived for two identical rigid spherical surfaces [l9| . 
The theory relates the pressure P generated in the liquid 
to the relative displacement of the two particles as: 



d , o, . dPir\). „ dh 
— (r x H 3 n — ^ = 12m — , 
dr\ dri dt 



(2) 



where, H(n) = ho + 2rf/d is the distance between the 
two surfaces at radial distance r\ from the center [see 
Fig. H(a)]. If the particles are completely immersed in 
the fluid, the expression for the viscous force is: 



3 2 1 dh 
'8^ h dt ' 



(3) 



For a liquid bridge, there is a correction coefficient JT^] 
to Eq. (3) and the expression for the force is: 



F V i s = --TTvd 2 (l - 



ho 



H{R) ho dt 



(4) 



where, R is the radius of the contact area. This equa- 
tion can be related to the volume of the liquid bridge 
using V — J Q R 2nriH(ri)dri. Pitois et. al. experi- 
mentally investigated the effect of viscosity of the fluid 
on the properties of a liquid bridge between two moving 



To study the segregation of the particles in the pile, 
we use different colors for the two kinds of particles. 
Fig. ^|(a) shows a pile after dry bidisperse granular mate- 
rial has been poured inside the silo. Here the small par- 
ticles appear white and the large particles appear black. 
Thus strong size separation is observed as the granular 
matter flow down the inclined surface. Two main mecha- 
nisms are important in determining the observed spatial 
distribution: (i) the void filling mechanism where smaller 
particles percolate through the larger particles and are 
thus found at the bottom of the flow [O, and (ii) the 
capture mechanism where the smaller particles which are 
more sensitive to surface fluctuations are stopped at the 
top of the pile before the larger particles |2(J . Segrega- 
tion is observed to vanish when a small amount of fluid is 
added. Fig. ^(b) shows the situation where Vf = 6x 10~ 3 , 
which corresponds to less than 1% by volume of water in 
the system. Thus a very small amount of fluid is observed 
to prevent segregation from occuring. 

To parameterize the extent of segregation, a histogram 
of the ratios of the two types of particles in a 3.5 mm 2 
area is measured using the light intensity. The light in- 
tensity is a monotonic function of the density ratio of the 
particles and is measured by using pre-determined weight 
ratios of particles in a separate series of calibration exper- 
iments. The histogram of local particle density ratio P w 
of the two kinds of particles in the pile is thus obtained 
and is plotted in Fig. ||(c),(d). The peaks in distribution 
can be fitted to a Gaussian and the mean value is used 
to determine the most common particle ratios a and b. 



The segregation parameter is defined as s 



(b-a) 
100 



which 



has a value between and 1 and describes the extent of 
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segregation. 



A. Dependence of the segregation on the size ratio 

s for three different bidisperse granular mixtures BD- 
1, BD-2, BD-3 (see Table |) are plotted in Fig. 0(a) as 
a function of Vf of water, s depends on the size ratio 
r of particles, and goes to zero at lower Vf for smaller 
r. However, segregation can persist even in the presence 
of fluid although the extent is considerably smaller. A 
phase diagram of the observed segregation as a function 
of r and Vf can be found in Ref. |14[| . 



B. Discussion of segregation on r size ratio 

The observed properties can be described in terms of 
the differences in the strength of interactions between 
particles due to difference in the sizes. In the presence of 
a wetting fluid, there are three possible interaction be- 
tween the particles: interactions between the (i) small- 
small particles S— S, (ii) small-large particles S — L, and 
(iii) large-large particles L — L. In addition, the particles 
and clumps will also interact with the inclined surface on 
which they flow. 

At small Vf the fluid is observed to selectively coat the 
small particles forming clumps of small particles due to 
the cohesive forces discussed earlier. Most of these inter- 
actions are in form of S — S. There is less percolation 
of small particles through larger particles thus reducing 
s. By increasing the amount of fluid further, the size of 
clumps increases, the interaction of S — L also become 
important, resulting in even smaller s. Finally, above a 
certain Vf ~ V c all types of interaction can be present, 
clumps are formed with both kinds of particles and seg- 
regation saturates. For large r clumps of S — L appear 
but are not strong enough and segregation decreases but 
saturates to non-zero values. 



C. Layering instability 

Fig. ^|(a) shows the intermediate situation where par- 
tial segregation is observed. We obtain the azimuthal 
correlation function g((j>) of the particle density inside 
the pile to characterize the layering where g(cj>) is given 
by: 

x - I(x,y) x I(x + r 2 cos((f>),y + r 2 sin ((/>)) 

gW = 2^ WW (5) 

x,y,r 2 y 

where, I(x,y) is the particle ratio at position (x,y) and 
T2 and <j> define the distance from the point where corre- 
lation is calculated. The correlation function was cal- 
culated for a square region in the middle of the pile. 
The peaks were observed to decrease in amplitude and 



are hidden by the noise if the whole pile was used for 
averaging. Fig. |l(](a) shows g(4>) for the image shown 
in Fig. |](a). The correlation function has two distinct 
peaks corresponding to the angles parallel to the surface. 
The peaks are observed to decrease in amplitude as lay- 
ering fades with increased Vf. Finally for the complete 
mixture we observe the flat line shown in Fig. |l^(b) for 
Vf = 6 x 1CP 3 . The correlation function does not show 
peaks for larger size ratios consistent with the visual ob- 
servation. 

Although the layering is not as periodic as in dry mix- 
tures of rough and smooth particles some stratifi- 
cation is clearly observed below V c . We believe that the 
layering is related to the selective coating and clump- 
ing of smaller particles (see Fig. |9|). A clump of small 
particles effectively behaves as a particle with a rough 
surface which can be greater than the individual large 
particle. Makse et al. have shown that stratification 
can occur when two species with different surface rough- 
ness are present in granular flows when the size of the 
rough particles exceeds the smoother particles Jl2j . The 
features observed in Fig. ^ appear to be related to this 
mechanism, although the increased stick-slip nature of 
the flow at the surface brought about by the addition 
of the fluid makes the layering aperiodic. The strongest 
layering is observed for BD-2 and decreases at higher r 
where the clumps of small particles stay smaller than the 
larger beads. At higher Vf, the average size of clumps in- 
creases and the interaction between small and large par- 
ticles also become important. Thus clumps with small 
and large particles and layering disappears. 



D. Viscosity dependence of the segregation 

Next we examine the effect of the viscosity of the fluid 
on the extent of segregation. As in the measurements for 
r , we change the viscosity of the fluid by using mixtures 
of water and glycerol to keep the surface tension almost 
constant, s corresponding to water and glycerol is plot- 
ted in Fig. ||. The segregation is observed to decrease to 
zero at a lower value of Vf in the case of higher glycerol, 
which has v a thousand times that of water. In fact, al- 
though some incipient amount of segregation is observed 
for large r in the case of water, segregation is observed 
to completely dissappear in the case of glycerol for all 
measured r. 

In Section III.C we discussed the strength of the vis- 
cous forces that may be important in determining perco- 
lation. For the particles [d ~ 1.2 mm) used in the exper- 
iments, and for glycerol or water as interstitial fluid, the 
saturation force due to surface tension is of the order of 
1(T 4 N. Therefore the viscous force is greater than the 
capillary force in the case of glycerol. The viscous force 
damps velocity fluctuation as it increases with relative 
velocity between particles. Because velocity fluctuations 
and percolation are required for segregation, s is there- 
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fore observed to be lower at higher v for similar T. Since 
the linear increase of 9 r , and decrease of s occurs be- 
low a similar V c , and both of these quantities saturate 
above this volume fraction, the cohesion due to the fluid 
appears to have the same effect on both 9 r and s. 



V. EXPERIMENTS WHEN PARTICLES ARE 
IMMERSED IN THE FLUID 

In this section, we study the angle of repose and the 
progress of segregation when the bidisperse glass parti- 
cles are poured into a silo filled with various liquids with 
different v. The procedure and the system is similar to 
the previous sections except that in these experiments, 
we first fill the silo with the fluid and then pour the 
dry granular matter. In this case only viscous forces are 
present and capillary forces are not present because liq- 
uid bridges are absent. 

Because the terminal velocity for particles decreases 
with i/, they take a significantly longer time to drain in 
a viscous fluid compared to that in air and Q < 0.02 
g s _1 in glycerol. As v is increased, the particles are 
increasingly observed to deviate from a downward tra- 
jectory and are deflected further down the slope even 
before coming in contact with the surface. The reason 
for the deflection of the particles can be qualitatively un- 
derstood as follows. For a particle moving with velocity 
v, a boundary layer develops in the fluid which is pro- 
portional to yj ud/pv. This estimate is essentially from 
dimensional arguments . As the distance between the 
particle and the surface becomes comparable to the size 
of this boundary layer, the particles feel a net horizontal 
force similar to the viscous forces discussed earlier due to 
fluid expulsion. These forces tend to deflect the particles 
further down the slope. 

We observed 9 r is ~ 24° ± 1° for mono and bidisperse 
glass particles for not only water and glycerol but also 
using polybutene (see Table ||) . Unlike the case of par- 
tial Vf, no systematic dependence of 9 r is observed on 
the viscosity. Furthermore, the observed 9 r in the fluid 
is very similar to that observed in air. This appears to 
suggest that the viscous forces discussed in Section III.C 
may not be relevant in the absence of liquid bridges. Al- 
though it must be noted that the experiments conducted 
in the regime where the particles are mixed with a small 
volume fraction of fluid is somewhat different in nature 
than these experiments because of the presence of the 
boundary layer. 

Experimental work investigating liquid bridges be- 
tween moving spheres 17 2^£^] has shown that bridge 
rupture distance and time between particles increases 
with viscosity. It is therefore possible that the increases 
in 9 r with viscosity in the case of partial Vf is related 
to the increase in the number and change in structure of 
liquid bridges with viscosity. 

Fig. 11(a) shows an image of the granular pile after 
bidisperse particles (BD-2) were poured into the silo filled 



with glycerol. Separate series of experiments to calibrate 
the intensity of light scattered by particles to the num- 
ber of particles were first conducted when particles were 
completely immersed in the various fluids. Using the 
same procedure as before, we measured the histogram of 
light intensity to determine s [see Fig. [n](b)]. Comparing 
Fig. |ll](a) to Fig. ^(a) shows a significant drop in s, when 
particles are poured in glycerol instead of air. The s is 
plotted in Fig. [T^(a) as a function of v for three different 
r. The segregation is observed to decrease as a function 
of v and drops to zero at high v which depends on r. 

The observed decreases in segregation with v can be ex- 
plained as follows. The boundary layer washes away the 
details of the surface roughness of the pile at higher v. 
Thus the capture mechanism which is sensitive to surface 
roughness decreases in importance with v. Furthermore, 
velocity fluctuations at higher v are damped out as parti- 
cles reach terminal velocity over a short distance leading 
to a reduction in percolation of particles There- 
fore the segregation decreases because the two mecha- 
nism which cause segregation in dry granular matter di- 
minish in strength. 

It is well known that sedimentation of different sized 
particles (for example in lakes) cause layering of different 
size particles to occur at the bottom because the terminal 
velocity reached by a spherical particle is proportional to 
\fd. This does not occur in our geometry because parti- 
cles are constantly fed into the system from the reservoir. 



VI. CONCLUSIONS 

In conclusion, the role of the size ratio of the parti- 
cles, volume fraction, surface tension and viscosity of the 
fluid on the extent of segregation and the angle of repose 
of a granular pile is clarified based on experiments and 
physical arguments. The angle of repose of the granu- 
lar matter is observed to increase sharply as the volume 
fraction of the fluid is increased and then saturates. The 
saturation occurs at a higher Vf than estimated by us- 
ing only particle roughness ideas. We observe that the 
viscosity of the fluid has a significant effect on the angle 
of repose of the pile. A sharp reduction of segregation is 
observed in the granular flow when a small volume frac- 
tion of fluid is added. The sharp changes in the angle of 
repose and segregation occurs over similar volume frac- 
tions suggesting that cohesivity has the same effect on 
both properties. The experiments point to a need for 
examining the role of the number of liquid bridges be- 
tween particles in determining the angle of repose. Our 
experiments appear to indicate that the changes in the 
number of liquid bridges and their structure with volume 
fraction, surface tension and viscosity of the fluid may 
be very important in determining the properties of wet 
granular matter. 
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and G. McKinley for bringing Ref. to our attention. 
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FIG. 1. Schematic diagram of the experimental setup. The 
reservoir containing the initial mixture is pulled up using a 
stepper motor in order to keep the height through which the 
mixture falls constant. 



FIG. 2. Images of a pile of dry and wet monodisperse par- 
ticles, (a) V f = 0, d = 0.5 mm, (b) V) = 10~ 2 , d = 0.5 mm, 
(c) V f = 0, d = 3.1 mm, (d) Vf = 10~ 2 , d = 3.1 mm. r 
is observed to be greater for the wet case and the surface of 
the pile of small particles is observed to be rough compared 
to the larger particles. 
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FIG. 3. (a) Size dependence of the r for d — 0.5 mm, 
d — 1.2 mm, and d — 3.1 mm glass beads. Note that 6 r is 
lower for larger d. (b) The saturation value of the 9 r of bim- 
ixtures are in between the 6 r of their components, (c) The 
mean square of the deviation of the surface from a straight 
line. The smaller particles are observed to stick together more 
readily, and make more clumps giving rise to a rougher surface 
and higher \ ■ 
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FIG. 4. (a) r for monodisperse particles as a function of 
Vf, for water {y — 10~ 3 N s m~ 2 ), glycerol (y = 1.5 N s 
m~ 2 ) and water-glycerol mixture (y — 0.12 N s m -2 ). (b) 6 r 
as a function of v for different bimixture and monodisperse 
particles for V f — 24 x 1CT 3 . 



FIG. 5. (a) Schematic of a liquid bridge between two par- 
ticles, (b) Image of a liquid bridge between two 0.1 mm par- 
ticles corresponding to Vf = 24 x 10~ 3 . 



FIG. 7. The extent of segregation s observed in bidisperse 
particles as a function of Vf of water for various size ratios r. 
Segregation decreases and saturates, when Vf is greater than 
V c . Segregation can persist if r is large. 
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FIG. 6. Image of the granular pile after bidisperse glass 
beads are poured into the silo (r = 2.4). The small particles 
are white and the large particles are black, (a) Vf = 0, (b) 
Vf = 6 x 10~ 3 . Strong segregation is observed for the dry 
case, the segregation is drastically reduced when less than 1% 
of water is present in the mixture, (c)-(d) The histogram of 
local particle density ratio of the two kinds of particles in the 
pile. The peaks in the distribution can be fitted to Gaus- 
sian and the mean value used to determine the most common 
particle ratios. 



FIG. 8. Viscosity dependence of the s and 8 r for r = 6.2 as 
a function of Vf . Note that the saturation value V c is similar 
for s and 9. 
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FIG. 9. (a) Strong layering is observed at Vf below V c for 
r = 2.4 (here Vf = 1.2 x 10~ 3 ). (b) The fluid is observed to 
selectively coat the smaller particles, resulting in clumps of 
small particles. In this case the clumps of small particles are 
larger than the larger particles. 
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FIG. 10. Azimuthal correlation function of the density of 
the particles for (a) V) = 1.2 x 1CT 3 and (b) V) = 6 x 1CT 3 . 
The correlation function show two distinct peaks at the an- 
gles parallel to the surface. The peaks fade for higher Vf . The 
strongest layering is observed for r — 2.4. 



FIG. 12. (a) The extent of segregation s(v) is observed to 
decrease as a function of v of the fluid, (b) 8 r as a function v 
is approximately constant. This is in contrast to observations 
for partial Vf where B r is observed to increase. 



Monodispcrsc rf(mm) 

0.1 ± 0.1 
0.5 ± 0.1 
0.9 ± 0.1 
1.2 ± 0.1 
3.1 ± 0.1 



mono-1 
mono- 2 
mono-3 
mono-4 
mono- 5 



23.5° ±0.5° 
24.0° ±0.5° 
24.0° ±0.5° 
24.0° ±0.5° 
25.5° ±0.5° 



Bidisperse 


d large 


d small 


r 


9 r 


BD-1 


mono-4 


mono-3 


1.3 


23° ±0.5° 


BD-2 


mono-4 


mono- 2 


2.4 


23° ± 0.5° 


BD-3 


mono- 5 


mono- 2 


6.2 


25° ±0.5° 


BD-4 


mono-5 


mono-1 


31 


24° ±0.5° 



TABLE I. Monodisperse and bidisperse mixtures of spher- 
ical glass particles used in the experiment, d is the diameter 
of particles, r the size ratio, and 9 r the angle of repose of the 
dry particles. 



FIG. 11. Image of the granular pile after bidisperse glass 
beads (r = 2.4) are poured into a silo filled with glycerol. 
Compare the extent of segregation to Fig. |^(a), where air 
is the interstitial medium. The segregation is drastically re- 
duced when particles are poured into a more viscous medium, 
(b) The histogram of local particle density ratio of the two 
kinds of particles in the pile. 
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-4 -2 n 2 

10 10 10 10 

V (N s m ) 



Fluid v 


(Nsm~ 2 ) 


p (kgm A ) 


T (Nnr 1 ) 


Water 


0.0010 


997 


0.07 ±0.003 


Water-glycerol- 1 


0.0098 


1076 


0.07 ±0.003 


Water-glycerol-2 


0.1190 


1108 


0.07 ±0.003 


Glycerol 


1.5 


1126 


0.07 ±0.003 


Polybutene L-50 


0.22 


844 


0.03 ±0.003 


Polybutene H-300 


30 


892 




Polybutene H-300 


100 


892 





TABLE II. The fluids used in the experiments, v is the 
viscosity, p the density and F is the surface tension of the 
fluid. All data corresponds to 25° C except in case of poly- 
butene H-300, where the lower v is obtained by heating the 
tank to 50° C. Water-glycerol- 1 and water-glycerol-2 contain 
60% and 88% glycerol by weight, respectively. 
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